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ABSTRACT 


In  experiments  going  back  to  the  first  free  electron  laser  (FEL)  oscillator  at  Stanford,  the 
measured  width  of  the  desynchronism  curve  is  often  significantly  greater  than  predicted  by 
theory  and  two-dimensional  (2D)  simulations  in  (z,  t).  The  results  of  new  four-dimensional 
(4D)  simulations  in  (jc,y,z,f)  show  that  this  difference  can  be  explained  by  the  effects  of 
diffraction.  When  the  light  is  artificially  constrained  to  remain  in  the  cavity  fundamen¬ 
tal  mode,  2D  and  4D  simulations  give  similar  results,  but  when  the  light  is  allowed  to 
self-consistently  develop  higher-order  modes,  the  4D  simulations  give  different  results  that 
agree  better  with  experiments.  The  results  of  new  4D  simulations  also  show  the  effects 
of  emittance  versus  electron  beam  energy  and  mirror  shift  versus  mirror  tilt  on  extraction. 
Analysis  of  these  results  examine  the  robustness  of  FEL  designs. 
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CHAPTER  1: 
Introduction 


1.1  A  Brief  History 

The  first  free  electron  laser  (FEL)  was  built  at  Stanford  University  and  lased  in  1976  [1]. 
This  FEL  was  initially  proposed  by  Madey  [2]  in  1971  and  was  based  on  initial  theory 
and  experiments  by  Motz  [3],  [4].  Stanford  then  went  on  to  build  and  operate  the  first 
oscillator  FEL  in  1977  [5].  Since  1976,  FELs  have  been  established  across  the  world  [6]- 
[8]  with  varying  wavelengths  and  power  levels.  The  most  powerful  FEL  was  operated 
at  Jefferson  Lab  in  2001  [9]  with  an  average  power  of  14kW  as  part  of  a  United  States 
Navy  (USN)  program  to  design  and  build  a  MW  class  FEL  for  ship  defense.  The  National 
Academy  of  Sciences  [10]  created  a  report  in  2009  detailing  the  issues  involved  in  the 
design  of  a  MW  class  FEL.  The  Office  of  Naval  Research  (ONR)  funded  an  Innovative 
Naval  Prototype  (INP)  in  2010  with  the  goal  of  designing  a  prototype  for  a  lOOkW  FEL. 
Boeing  was  selected  for  this  project;  however,  it  was  canceled  in  2012  as  the  Navy  decided 
to  focus  on  solid  state  lasers  (SSL)  for  their  near  term  deployment  capabilities. 

For  Naval  applications,  directed  energy  (DE)  weapons  represent  an  opportunity  to  gain 
a  significant  advantage  in  both  offensive  and  defensive  capabilities.  A  SSL  has  already 
been  deployed  on  a  surface  combatant;  however,  SSLs  are  limited  in  power  and  to  a  spe¬ 
cific  wavelength.  Previous  research  at  Naval  Postgraduate  School  (NPS)  [11]  has  shown 
that  there  are  wavelengths  outside  of  the  typical  SSL  operation  that  enhance  propagation 
through  the  atmosphere  in  certain  maritime  environments.  FELs  are  capable  of  lasing  at 
these  wavelengths  with  better  beam  quality  and  at  potentially  higher  powers  that  provide 
better  utility  in  their  application.  Continuing  research  regarding  FELs  is  critical  for  the 
USN  to  maintain  an  advantage  in  DE  weapons. 


1.2  Overview  of  Thesis  Research 

The  DE  group  at  NPS  has  been  developing  FEL  modeling  programs  for  decades.  Recently, 
it  developed  a  four-dimensional  model  that  follows  the  evolution  of  the  electron  and  opti- 
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cal  pulses  in  the  transverse  and  longitudinal  directions.  Previous  models  have  been  used  to 
simulate  FELs,  but  some  discrepancies  with  comparisons  to  theoretical  and  experimental 
results  existed.  This  thesis  work  looks  at  these  discrepancies  using  the  new  model  for  the 
original  Stanford  FEL  and  the  Jefferson  Laboratory  14kW  FEL.  Research  into  the  design 
of  high  power  FELs  has  also  been  conducted  by  the  DE  group,  but  this  was  using  previ¬ 
ous,  less  sophisticated  models.  This  thesis  expands  this  work  by  using  the  new  model  to 
determine  the  robustness  of  FEL  designs. 

1.3  Outline 

Chapter  2  derives  the  equations  that  constitute  the  base  of  FEL  theory  as  well  as  discusses 
in  detail  the  concepts  of  pulse  lethargy  and  desynchronism  [12],  [13].  Chapter  3  discusses 
the  codes  used  to  model  FELs  and  the  output  from  the  simulations.  Chapter  4  establishes 
the  design  parameters  of  a  high  power  FEL.  Chapter  5  lays  out  the  results  of  the  simulations 
and  Chapter  6  discusses  the  results.  Chapter  7  provides  a  summary  of  the  work  as  well  as 
a  discussion  regarding  future  research. 
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CHAPTER  2: 

Theory  of  Free  Electron  Lasers 


2.1  Component  Description  of  an  Energy-Recovery  Linac 
Free  Electron  Laser 

The  components  of  an  energy-recovery  linear  accelerator  (linac)  FEL  can  be  seen  in  Figure 
2.1.  The  first  component  of  an  energy-recovery  linac  FEL  is  the  electron  injector.  The 
injector  can  produce  the  electrons  by  different  methods  (e.g.,  by  a  thermionic  process  or 
by  photoemission).  These  electrons  are  produced  in  a  pulse  structure  with  a  typical  bunch 
charge  (q^)  of  0.1  nC  and  duration  (//,)  of  1  ps.  A  typical  value  of  the  electrons’  kinetic 
energy  after  the  injector  is  5  MeV.  The  electrons  are  then  directed  into  the  linac,  which 
often  consists  of  radio  frequency  (RF)  powered  accelerators.  For  such  accelerators,  an  RF 
signal  is  applied  to  the  accelerator  cavities  to  generate  an  alternating  electric  field  inside 
them.  The  alternating  electric  field  continually  accelerates  the  electrons  throughout  the 
cavities  as  the  electrons  “surf”  the  RF  wave.  The  electrons  in  the  RF  linac  are  accelerated 
to  a  typical  energy  of  100  MeV. 


resonator  mirror 


Figure  2.1:  Diagram  of  the  components  of  an  Energy  Recovery  Linac  FEL 
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The  electrons  are  transported  between  each  component  via  an  electron  beam  line.  This  line 
is  an  evacuated  pipe  that  has  a  radius  of  ~  1  cm  and  is  surrounded  by  various  magnets 
and  detectors.  Solenoid  and  quadrupole  magnets  are  used  to  focus  the  beam  so  that  the 
electrons  stay  collimated.  Dipole  magnets  are  also  used  to  steer  the  beam  at  various  places. 

The  next  component  is  the  undulator.  It  is  comprised  of  permanent  magnets  or  electromag¬ 
nets  (with  peak  magnetic  fields  of  approximately  1  T)  arranged  in  such  a  way  to  produce  an 
alternating  magnetic  field.  As  the  electrons  move  through  the  undulator,  they  wiggle  back 
and  forth  due  to  the  Lorentz  force  and  emit  photons.  These  photons  are  often  stored  in  an 
optical  cavity  consisting  of  two  mirrors:  a  highly  reflecting  (~  100  %  reflective)  mirror  and 
a  partially  transparent  (~  1  —  10  %  transparent)  mirror.  The  reflected  optical  pulse  inside 
the  cavity  arrives  at  the  undulator  at  approximately  the  same  time  as  an  electron  pulse.  The 
electron  pulse  then  amplifies  the  optical  pulse  through  the  production  of  additional  photons. 
This  process  continues  over  many  passes  until  the  gain  from  the  electron  pulse  equal  the 
losses  (due  to  the  partially  transparent  mirror),  which  is  referred  to  as  saturation.  The  laser 
output  power  depends  on  the  saturated  power  inside  the  cavity  and  the  mirror  transmission. 

The  electrons  are  then  recirculated  back  to  the  RF  linac  using  dipole  magnets.  The  electrons 
are  now  injected  180°  out  of  phase  with  the  RF  fields  that  causes  the  electrons  to  decelerate 
back  to  near  their  original  injection  energy  of  ~  5  MeV.  This  takes  energy  away  from  the 
electrons  and  transfers  it  to  the  RF  field  (energy  recovery).  Once  the  electrons  exit  the  RF 
linac,  they  are  directed  to  a  beam  dump  where  they  impact  shielded,  conducting  material. 
This  method  of  recycling  the  electrons  greatly  enhances  the  overall  “wall  plug”  efficiency 
and  reduces  the  generated  radiation  at  the  beam  dump. 


2.2  FEL  Resonance  Condition 

Resonance  refers  to  the  condition  in  which  optimum  energy  exchange  occurs.  Specifically, 
Colson  [14]  describes  this  as  the  case  where  as  an  electron  passes  through  one  undulator 
period,  one  wavelength  of  light  passes  over  it,  which  is  important  so  that  the  electron  sees 
the  correct  electric  field  as  it  traverses  the  undulator.  As  the  electron  moves  through  the 
undulator,  there  are  two  speeds  that  must  be  considered:  the  speed  c  of  the  photon  and  the 
speed  /3 zc  of  the  electron  with 
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(2.1) 


Pz  = 


c 


(1  +K2) 
If  ’ 


where  K  =  eBAo / Inmc2  is  the  undulator  parameter,  7=1/ ^/l  —  (v2/c2),  <?  is  the  charge  of 
the  electron,  B  is  the  root  mean  square  field  strength  over  each  period,  Aq  is  the  undulator 
period,  and  m  is  the  mass  of  the  electron.  The  approximation  made  above  will  be  explained 
in  the  next  section.  The  time  for  an  electron  to  travel  one  undulator  period  is  At  —  Ao//3 zc. 
In  that  time,  the  photon  will  move  ahead  of  the  electron  by  a  distance  of  (c  —  /Tc)  At,  which 
we  equate  to  the  optical  wavelength  as  a  condition  for  resonance  where 


A  =  (c  —  fizc)  At 


(1  ~  Pz)  Aq 

P z 


(2.2) 


Placing  (2.1)  into  (2.2)  gives  the  approximate  wavelength  the  FEL  will  lase  at: 


A0(l+£2) 

If 


(2.3) 


Equation  (2.3)  shows  one  of  the  advantages  of  an  FEL  over  a  solid  state  laser  since  the 
FEL  wavelength  is  tunable  by  manipulation  of  either  K  or  y.  The  wavelength  can  also  be 
manipulated  by  adjusting  Ao,  however  this  is  usually  not  done  as  it  is  easier  to  manipulate 
the  magnetic  field  strength  or  the  electron’s  energy  coming  out  of  the  accelerator. 


2.3  The  Pendulum  Equation 

An  electron  traveling  though  a  helical  undulator  along  the  z-axis  is  subjected  to  a  magnetic 
field  described  as: 


B  =  fi(cos(Aoz),sin(Aoz),0), 

where  B  is  the  magnetic  field  amplitude  and  ko  =  27t/Ao  is  the  undulator  wavenumber.  The 
magnetic  and  electric  fields  of  the  laser  are  given  by 
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fi  =  £(sin(vP),cos(vP),0), 


(2.4) 


E  =  £,(cos(vP),  —  sin(vP),0),  (2.5) 

where  E  is  the  electric  field  amplitude  (in  cgs  units),  =  kz  —  CQt  +  (j)  is  the  optical  phase, 
k  —  2n/X  is  the  wavenumber,  A  is  the  optical  wavelength,  and  co  is  the  optical  frequency. 
The  relativistic  Lorentz  force  equations  are 


djyp) 

dt 


e 

me 


(E  +  pxB), 


(2.6) 


(2.7) 


(2.8) 


where  v  =  [3  c  is  the  electron  velocity. 

The  motion  of  the  electron  is  determined  by  the  combination  of  the  undulator  and  laser 
fields.  Placing  (2.4)  and  (2.5)  into  (2.6)  gives: 

=  — ~—[E{\  — /3z)(cos(vP),  —  sin(vP))  +  /3;5(—  sin(fcoz),cos(£oz))],  (2-9) 
clt  me 

=  -  —  [E(f}xc> os(lP)  -/3vsin(vP))  +5(/3^sin(A0z)  -/3vcos(fc0z))],  (2.10) 

where  /3 j_c  is  the  transverse  velocity  of  the  electron.  Since  the  electrons  are  relativistic, 
/3-  «  1,  which  means  (2.9)  can  be  simplified  then  integrated  to 
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(2.11) 


Ai  = - ( cos(koz),sm(k0z),0 ), 

r 

where  we  assume  that  the  electron  is  injected  into  its  orbits  such  that  the  constants  of 
integration  are  zero. 

To  determine  the  electron’s  energy  exchange  with  the  optical  field,  we  place  (2.11)  into 
(2.7)  to  get 


7~~t~  =  — — E[&cos('P)-j3ysin('F)]  =  ^-cos  (£  +  0),  (2.12) 

at  me  ■  Ymc 

where  the  electron  phase  is  £  =  (k  +  ko)z  —  (Ot.  The  variables  ko,  k,  and  co  are  constant 
which  implies  that  £  describes  the  electron’s  position  on  the  scale  of  an  optical  wavelength 
since  £  kz(t).  When  cos(£  +  <j>)  >0,  the  electron  energy  increases  and  y  >  0;  but  if 
cos(£  +  (j>)  <0,  the  electron  energy  decreases  and  y  <  0  (the  desired  state  as  it  leads  to 
optical  gain).  Sufficient  energy  exchange  between  the  electron  and  laser  beam  requires 
that  cos(£  +  0)  not  oscillate  quickly  over  time  L/c/ 3Z,  where  L  =  NX o  is  the  length  of  the 
undulator  and  N  is  the  number  of  undulator  periods. 

Now,  we  can  define  the  electron  phase  velocity  as 

v  =  —  =  L[(k  +  ko)  pz  —  k]  =  £, 
c 

where  (...)  =  indicates  a  derivative  taken  with  respect  to  dimensionless  time  x  =  ct/L. 
The  parameter  r  goes  from  0  to  1  as  the  electron  traverses  the  undulator. 

Notice  that  if  v  =  0,  j8z  =  kc/ (k  +  ko),  which  is  equivalent  to  the  resonance  condition  (2.2). 
Thus  v  =  0  corresponds  to  resonance,  where  the  energy  exchange  between  electrons  and 
light  is  maximized.  Substituting  (2.11)  into  (2.8)  gives 
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1- 


Pz  = 


l+^y 
r2  ) 


Using  the  relativistic  limit  where  y  1,  this  can  be  expanded  as 


Pz 


1  +K2 
2  y2 


Throughout  the  electron  beam,  the  initial  electron  phase,  £(0)  =  Co  =  (k  +  ko)zo  ~  feo> 
depends  on  the  initial  electron  position  zo  and  is  therefore  different  for  each  electron.  An 
inspection  of  (2.12)  shows  that  approximately  half  of  the  electrons  give  energy  to  the  laser 
beam  and  the  other  half  take  energy  from  the  laser  beam.  This  phenomenon,  with  half  the 
electrons  moving  faster  and  the  other  half  moving  slower,  causes  the  electrons  to  bunch  and 
radiate  coherently  due  to  the  fact  that  it  occurs  on  the  scale  of  A . 

If  we  consider  v  =  L[(k  +  ko)Pz]  where  j8z  «  1  —  (1  +K2)/ 2y2,  and  using  the  resonance 
condition  (2.3)  with  the  assumption  y  1,  then  we  find  that 


o  4kNy 

v  = - - 

y 


(2.13) 


Placing  (2.12)  into  (2.13)  gives 


v  —  ^  —  \ci\cos(^  +  (j)),  (2.14) 

where  \a\  =  AnNeKLE /y2mc2  is  the  dimensionless  laser  field  amplitude.  Equation  (2.14) 
has  the  form  of  the  classical  pendulum  equation  and  describes  how  the  electrons  evolve 
in  phase  space  (£,v).  An  inspection  of  (2.14)  shows  that  when  \a\  S>  K  the  fields  of  the 
laser  are  strong,  which  causes  the  electrons  to  bunch  together;  however,  if  \a\  <U  K  then  the 
fields  of  the  laser  are  weak  and  do  not  cause  significant  bunching.  The  electrons  follow 
FEL  phase  space  paths  described  by 
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v2  =  Vo  +2\a\ [sin(C  +  0)  -  sin(£0  +  0)], 


where  £o  and  Vo  are  the  initial  coordinates  of  the  electron  in  phase  space.  Figure  2.2  shows 
an  example  of  electron  evolution  in  phase  space  through  an  undulator.  Here,  you  can  see 
the  electrons  are  introduced  at  the  beginning  of  the  undulator  with  an  initial  v.  As  the 
electrons  progress  through  the  undulator,  they  begin  to  evolve  in  phase  space.  At  the  end 
of  the  undulator,  you  can  see  that  the  electrons  have  evolved  significantly  from  their  initial 
v  values  and  are  now  bunched  on  the  scale  of  an  optical  wavelength. 


electron  beam 


Figure  2.2:  Evolution  of  electrons  through  the  undulator  in  phase  space. 


2.4  The  Wave  Equation 

Maxwell’s  wave  equation  is 


=  <2'15) 

where  A(x,t)  is  the  optical  vector  potential  and  J±  is  the  perpendicular  component  of  a 
current  source.  The  electric  and  magnetic  fields  are  determined  using 
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1  dA 

E  = - ~~T7 i 

c  at 

B  =  VxA. 

For  a  laser,  we  can  assume  the  electric  field  envelope  is  slowly  varying  with  distance  on 
the  scale  of  an  optical  wavelength,  A  —  2 n/k,  and  in  time  when  compared  to  the  optical 
period,  2 n/ (O  —  X/c  [14].  Therefore,  this  motivates  us  to  write  the  optical  vector  potential 
as 


A(x,t) 


E(x,t) 

k 


£e 


ia 


5 


where  a  =  kz  —  cot  is  the  phase  of  the  carrier  wave,  E  =  \E\el *  is  the  complex  laser  field 
envelope,  and  e  is  the  polarization  vector.  Both  the  amplitude,  \E(x,t)\,  and  the  phase, 
<j)(x,t)  of  the  wave  envelope  are  slowly  varying,  therefore  we  can  rewrite  (14)  as 


£e 


ia 


k 


Vl  +  2ik 


E 


(2.16) 


where  V2_  =  +  d2 . 

Using  the  process  of  method  of  characteristics,  a  coordinate  substitution  can  be  made  such 
that 


E(z,t)  =  E'(u,t), 

where  u  —  z  —  ct.  Now,  we  can  write 

dE  _  dE'  du  _  dE' 
dz  du  dz  du  ’ 


(2.17) 
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since  du/dz  =  1.  Likewise, 


BE  _  BE'  du  dE'_  BE'  BE' 

Bt  Bu  Bt  Bt  u  Bt 

Using  (2.17)  and  (2.18),  we  can  write 


(2.18) 


BE  l  BE  _  1  BE' 

Bz  c  Bt  c  Bt 

Multiplying  both  sides  of  (2.16)  by  ke~m£*,  using  the  coordinate  substitution,  and  drop¬ 
ping  the  prime  notation  now  gives  the  paraxial  wave  equation 


V\+2ik 


E  = 


47 ik 


y,  -£*e~m. 


(2.19) 


The  source  current  can  be  written  as 


J  =  -ecJ2^5(3)(T-r/(t)), 

i 


(2.20) 


where  r,(/)  is  the  ith  electron  position  at  time  t.  The  contribution  to  the  transverse  current, 
from  the  transverse  motion  as  described  in  (2.11),  is 

j8 _l  =  -  —  (cos(k0z),sin(k0z),0)  =  Re  ( -—i£e~lk°z)  ,  (2.21) 

7  V  7  / 

where  £  —  4^(— /.  1 .0)  for  circular  polarization  from  a  helical  undulator.  Substituting 
(2.21)  into  (2.20)  and  then  (2.20)  into  (2.19)  gives 


2  ik  B 
c  Bt 


E 


—4  TCieKkp(x,t) 


(2.22) 
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where  p(x,t )  is  the  local  electron  density  and  (e_!l=/y)  represents  an  average  over  sample 
electrons  in  a  volume  element  dV  that  is  significantly  smaller  than  either  the  optical  or 
electron  pulse,  and 


f  ^8^3\x-^ri(t))dV  =  p(x,t)dV 

■'  i 

represents  the  number  of  electrons  in  the  volume  element. 

Here,  again,  we  will  use  the  dimensionless  time  parameter  T  =  ct/L  to  rewrite  (2.22)  as 


2  ik  d 
L  dz 


E 


— AnieKkp(x ,  t) 


(2.23) 


Multiplying  (2.23)  by  —AnN eKL2 / yfimc2k  now  gives  the  wave  equation  as 


a(x,z)  =  -j(e  ^), 

where  a  =  \a\et(^  —  AnNeKLE / y^mc2  represents  the  dimensionless  laser  field  amplitude 
and  the  dimensionless  FEL  current  is  represented  by  j  =  %n2Ne2K2L2p  jy\mc2  using  the 
assumption  that  y  ~  7o.  Defining.?  —  x( k/ 2L)1/2  andy  —y(k/2L)1^2  allows  us  to  now  write 
the  wave  equation  in  completely  dimensionless  form  as 


iL  -l7  d 

— V  i  H - 

2k  d  T 


a  {x,  t) 


where  for  simplicity  the  tildes  are  dropped  and  the  coordinates  (x,y,z,  r)  are  dimensionless. 

The  first  term  in  the  above  equation  describes  diffraction  of  the  optical  field.  If  we  consider 
the  case  where  the  electron  and  optical  beams  exactly  overlap  and  the  effect  of  diffraction 
is  small,  the  wave  equation  can  then  be  simplified  as 
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a  =  -j{e 


(2.24) 


where  the  coupling  between  the  laser  and  electron  beam  is  determined  by  j  and  the  electron 
bunching  by  (e~!<=).  For  small  amounts  of  bunching,  the  coupling  is  small.  When  the  value 
of  j  <  K.  the  coupling  is  small  but  when  j  3>  7T,  the  coupling  can  be  large. 

Using  the  fact  that  a  =  \a\elv ,  the  left  side  of  (2.24)  becomes 


a  =  \a\el ^  +  \a\^>iel<^. 


(2.25) 


Now,  placing  (2.25)  into  (2.24)  and  grouping  real  and  imaginary  terms  gives 


a\  =  -j(  cos(C +  <|>)), 


(2.26) 


^  =  j(sin(C  +  0)) 
\a\ 


(2.27) 


Whether  the  coupling  drives  the  optical  field  amplitude  \a\  or  phase  (f>  depends  on  the 
relative  values  of  £  and  (j),  as  can  be  seen  in  (2.26)  and  (2.27).  When  electrons  bunch, 
such  that  £  4-  (f>  — >  K,  the  value  of  (2.26)  increases.  This  increase  in  the  dimensionless  laser 
field  amplitude  in  turn  causes  the  electrons  to  bunch  even  more  via  the  coupling  described 
in  (2.14).  As  more  bunching  occurs,  this  causes  £  and  (j)  to  become  more  coupled,  thus 
increasing  the  value  of  (2.26).  This  process  is  a  kind  of  “feedback”  loop  between  the 
pendulum  and  wave  equations  for  the  FEL  [14].  This  loop  causes  exponential  growth 
initially,  until  saturation  occurs.  As  the  value  of  \a\  increases,  the  evolution  of  the  electrons 
in  phase  space  increases.  Once  the  value  increases  sufficiently,  £  +  0  — y  0,  which  causes  the 
electron  energy  to  increase,  therefore  taking  energy  from  the  laser  beam.  This  saturation 
effect  prevents  the  runaway  of  the  exponential  grow  in  gain. 
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2.5  Pulse  Lethargy  and  the  Detuning  Curve 

Following  the  phase  space  evolution  of  the  electrons  through  the  undulator,  as  described 
by  (2.14),  the  electrons  enter  the  undulator  unbunched  at  T  =  0.  As  the  electrons  progress 
through  the  undulator,  they  begin  to  bunch  until  they  reach  the  end  of  the  undulator,  T  =  1, 
where  the  electrons  have  bunched  such  that  they  significantly  amplify  the  optical  pulse. 
Recall  from  the  previous  discussion  on  the  resonance  condition  that  the  photons  “slip" 
ahead  of  the  electrons  by  a  distance  of  NX  over  the  length  of  the  undulator.  This  is  referred 
to  as  the  slippage  distance  [12].  This  slippage  causes  the  electrons  to  amplify  the  trailing 
end  of  the  optical  pulse,  causing  the  centroid  to  move  slower  than  c  as  shown  in  Figure  2.3. 
This  effect  is  called  lethargy  [13]  and,  over  many  passes,  will  cause  the  optical  pulse  to  fall 
behind  the  electron  pulse  and  the  amplitude  to  decay. 


Figure  2.3:  Depiction  of  pulse  lethargy  with  the  electron  pulse  in  red  and  the  optical  pulse  in 
blue. 


To  counteract  this  lethargy  effect,  the  spacing  between  the  mirrors  of  the  optical  cavity 
must  be  adjusted  to  shorten  the  optical  cavity.  This  process  is  called  desynchronization 
and  allows  the  photon  pulse  to  be  advanced  after  each  pass  [12].  If  the  cavity  length  is 
shortened  by  moving  one  mirror  a  distance  AS,  then  the  amount  of  desynchronize  can  be 
characterized  by  the  dimensionless  parameter 


d 


-2AS 

NX 


(2.28) 


The  plot  of  the  electron  extraction  rj  vs.  d  is  called  the  detuning  curve,  where  the  extraction 
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is  the  output  optical  power  over  the  input  electron  beam  power;  Figure  2.4  gives  an  example 
of  one.  From  this  example,  there  is  clearly  a  peak  that  will  provide  the  most  output  for 
the  laser,  however  this  region  can  become  unstable  in  the  event  of  a  slight  change  in  the 
mirror  separation.  Just  beyond  the  peak  is  a  region  that  provides  much  more  stability  with 
relatively  high  output.  At  large  values  of  d,  the  extraction  goes  to  zero  because  the  optical 
pulse  advances  too  far  ahead  of  the  electron  pulse  on  each  pass  [12]. 
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CHAPTER  3: 

FEL  Simulation  Codes  Overview 


The  modeling  programs  used  here  have  been  developed  by  professors  J.  Blau  and  W.B. 
Colson  at  NPS  [15].  All  the  programs  use  the  dimensionless  parameters  as  previously 
described.  The  two-dimensional  (2D)  program  follows  the  electron  and  optical  pulses  in 
the  longitudinal  (z)  direction  in  time  (t)  while  the  four-dimensional  (4D)  program  follows 
the  pulses  in  the  transverse  (x.y)  and  longitudinal  directions  in  time.  The  specific  4D  model 
used  is  for  short  pulses  (optical  pulse  length  comparable  to  the  slippage  distance  NX).  In 
the  model,  a  fourth-order  Runge-Kutta  method  is  used  to  evaluate  how  the  electrons  evolve 
in  phase  space  according  to  the  FEL  pendulum  equation  (2.14).  The  optical  field  propagates 
according  to  the  FEL  wave  equation  (2.24)  using  a  Fourier  transform. 

Figure  3. 1  is  an  example  of  the  first  portion  of  the  graphical  output  from  the  4D  model.  The 
shades  of  blue  represent  a  color  scale  with  the  darker  blue  indicating  a  lower  value  than  the 
lighter  blue.  The  bottom  row  of  plots  is  at  the  beginning  of  the  first  pass  and  the  top  row 
is  at  the  end  of  the  last  pass.  The  middle  row  of  plots  shows  the  evolution  over  n  passes 
through  the  undulator.  The  far  left  column  of  plots  is  labeled  a(x,  0,0)  and  represents  a 
cross-section  of  the  optical  field  amplitude  versus  x  where  y  —  z  —  0.  The  number  in  the 
top  right  comer  of  each  plot  gives  the  peak  value  of  \a\.  The  electron  beam  is  shown  in 
red  and  the  optical  beam  in  blue.  The  next  column  of  plots,  labeled  a(0,y,0),  is  nearly 
identical,  showing  a  cross-section  of  the  field  amplitude  versus  y  for  x  —  z  —  0.  The  third 
column  of  plots,  labeled  a(0,0,z),  shows  the  electron  pulse  and  optical  field  amplitude 
versus  the  longitudinal  direction  z,  at  x  =  y  —  0.  Initially,  the  electron  pulse  (red)  is  slightly 
ahead  of  the  optical  pulse  (blue),  but  over  the  length  of  the  undulator  the  electron  pulse  slips 
behind  the  optical  pulse  by  NX.  The  final  column,  labeled  P(0,0,  v),  shows  the  evolution 
of  the  optical  power  spectrum  at  x  =  y  =  0.  The  shift  in  the  centroid  of  the  optical  spectrum 
seen  here  illustrates  the  phenomenon  through  which  the  FEL  selects  the  best  wavelength  at 
which  to  operate,  as  it  evolves  from  weak  fields  to  saturation  [14]. 
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Figure  3.1:  First  sample  graphical  output  from  4D  model. 


The  second  set  of  plots  is  shown  in  Figure  3.2.  The  top  left  plot,  labeled  /(v,«),  shows 
the  evolution  of  the  electron  phase  velocity  distribution.  The  plot  next  to  that  shows  the 
final  distribution  of  electrons  in  phase  space  and  indicates  the  bunching  of  the  electrons. 
The  two  plots  in  the  upper  right  give  a  three-dimensional  (3D)  representation  of  the  optical 
wavefronts  at  the  left  mirror  (t  =  —0.1)  and  the  right  mirror  (t  =  1.1).  The  bottom  two 
plots  on  the  left  show  the  evolution  of  optical  power  P(n)  and  gain  G(n )  and  give  their  peak 
values  in  the  upper  right.  The  next  plot,  labeled  \c(m,p)\,  indicates  the  modal  composition 
of  the  optical  wavefront.  The  value  of  the  Hermite-Gaussian  coefficients  based  on  the 
cavity  modes  are  represented  by  a  blue  square  at  each  value  of  (m,  p)  using  the  previously 
mentioned  color  scale.  The  last  two  plots  give  a  2D  representation  of  the  optical  wavefront 
at  the  left  and  right  mirrors  with  the  peak  value  indicated  in  the  top  right. 
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Figure  3.2:  Second  sample  graphical  output  from  4D  model. 
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CHAPTER  4: 

Design  of  High  Power  FELs 


In  2013,  the  DE  group  at  NPS  published  an  article  [16]  discussing  the  design  requirements 
for  an  FEL  to  reach  an  average  power  of  1  MW.  The  DE  group  initially  started  with  a 
requirement  for  an  optical  wavelength  A  ~  1  /im  (chosen  for  atmospheric  propagation). 
Using  a  commercially  available  undulator  with  a  period  Ao  ~  3  cm  and  an  undulator  pa¬ 
rameter  K  ~  1  provides  the  required  wavelength  with  an  electron  beam  energy  E/,  ~  100 
MeV  (determined  using  (2.3)).  A  superconducting,  RF  (SRF)  linac  can  accelerate  elec¬ 
trons  with  an  accelerating  gradient  of  ~  10  MeV/m  which  gives  an  accelerator  length  of 
~10  m.  FELs  have  a  typical  extraction  77  ~  2%,  which  means  that  the  required  average 
electron  beam  power  is  Pf,  ~  50  MW.  Combining  that  with  the  electron  beam  energy  gives 
a  required  average  current  Iavg  ~  0.5  A;  this  can  be  achieved  with  a  bunch  charge  <://,  ~  1 
nC  at  a  pulse  repetition  frequency  of  500  MHz.  To  achieve  good  gain,  the  peak  current 
needs  to  be  Ipk  —  qb/th  ~  1  kA,  which  gives  a  bunch  length  th  ~  1  ps.  For  an  oscillator, 
the  extraction  is  approximately  r]  ~  1/2 N,  which  requires  N  ~  25  to  achieve  an  extraction 
r/  ~  2%,  producing  an  undulator  that  is  ~  0.75  m  long.  Most  of  the  design  parameters 
discussed  here  have  been  achieved  in  laboratory  FELs,  but  not  all  at  the  same  time  in  a 
single  FEL.  For  example,  a  few  electron  guns  have  produced  nC  bunch  charges,  but  not  at 
the  high  repetition  rate  listed  here,  so  their  average  current  is  much  less  than  the  goal  of 
~  0.5  A.  This  is  an  important  area  of  research  to  achieve  a  high-power  FEL. 

Previous  work  at  NPS  [17]  analyzes  the  robustness  of  both  an  oscillator  and  amplifier  FEL 
by  varying  emittance  (£x,y),  energy  spread  (A y/y),  mirror  shift  (A ym),  and  mirror  tilt  (A0m). 
The  Jefferson  Laboratory  14  kW  FEL  was  also  analyzed  for  comparison.  This  thesis  work 
focuses  on  the  oscillator  FEL  and  builds  upon  the  previous  work  at  NPS.  Normalized  emit¬ 
tance  is  defined  as  the  phase  space  area  of  the  electron  bunch  given  by  £XJ  =  yArxxAQxy 
where  A rXJ  is  the  rms  transverse  position  spread  of  the  electron  bunch  and  AQxy  is  the 
rms  angular  spread  of  the  electron  bunch,  as  shown  in  Figure  4.1  [18].  The  energy  spread 
refers  to  the  spread  in  energies  of  the  individual  electrons  in  the  electron  bunch.  A  typical 
accelerator  produces  an  energy  spread  Ay/ y  ~  0. 1%  — >■  1%. 
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The  mirror  shift  and  tilt  refer  to  misalignments  of  the  outcoupling  mirror,  which  could 
be  due  to  platform  vibrations  for  a  ship-based  FEL.  The  shift  is  simulated  by  moving  the 
mirror  up  (shifting  the  center  of  the  optical  beam  from  the  middle  of  the  mirror  toward  the 
bottom)  and  the  tilt  is  simulated  by  rotating  the  mirror  about  the  center  point,  which  can 
be  seen  in  Figures  4.2  and  4.3,  respectively.  The  mirror  shift  is  given  in  /im  and  the  tilt  is 
given  in  /irad.  The  mirror  shift  and  tilt  could  be  limited  to  these  magnitudes  using  active 
alignment. 


Figure  4.1:  Emittance  is  defined  by  the  area  of  the  electron  bunch  in  phase  space,  particle  angle 
(x')  versus  particle  position  (x),  from  [18] 
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Figure  4.2:  Mirror  Shift  (A ym) 
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Figure  4.3:  Mirror  Tilt  (A 6,„) 
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CHAPTER  5: 
Results 


5.1  Jefferson  Laboratory  FEL 


The  Jefferson  Laboratory  14  kW  FEL  was  modeled  using  both  the  2D  and  4D  programs  and 
the  specific  parameters  shown  in  Table  5.1.  This  FEL  was  chosen  to  be  modeled  because 
the  DE  group  at  NPS  has  previously  worked  with  the  team  at  Jefferson  Laboratory  and 
has  done  numerous  models  of  the  14  kW  FEL,  with  the  exception  of  using  the  recently 
developed  4D  program.  Many  experimental  results  for  this  FEL  have  been  measured  and 
published;  however,  a  desynchronism  curve  was  never  established.  This  does  not  allow  a 
comparison  of  the  modeling  results  with  the  experimental  results.  However,  an  examination 
of  the  modeling  results  does  provide  an  opportunity  to  visit  an  issue  discovered  by  Benson 
in  1985,  when  he  showed  that  the  width  of  the  desynchronism  curve  using  a  computer 
model  was  twice  the  length  of  the  experimental  desynchronize  curve  [19]. 

Figure  5.1  shows  the  desynchronism  curve  for  the  Jefferson  Laboratory  14  kW  FEL  using 
both  the  2D  and  4D  programs.  The  2D  results  are  in  blue  and  the  4D  results  are  in  red. 
The  experiment  achieved  a  peak  extraction  of  Jj  «  1.4%  which  compares  well  with  both  of 
the  curves,  but  the  2D  curve  is  approximately  half  the  width  of  the  4D  curve.  This  mirrors 
the  issue  described  by  Benson  where  the  4D  curve  presumably  represents  the  experimental 
data.  Closer  examination  of  the  4D  curve  shows  that  the  optical  wavefront  at  the  point  of 
peak  extraction  id  =  0.02)  is  96%  in  the  fundamental  cavity  mode.  Figure  5.2  shows  the 
electron  pulse  (in  red)  and  the  optical  field  amplitude  (in  blue)  versus  z  at  .v  =  y  =  0  and 
the  2D  representation  of  the  optical  wavefront  at  the  right  mirror  at  a  larger  desynchronism 
d  —  0.02.  An  examination  of  the  optical  wavefront  at  d  —  0.2  shows  that  it  is  now  only 
66%  in  the  fundamental  mode,  as  shown  in  Figure  5.3. 
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Electron  Beam: 

Energy 

Eb  =  115  MeV 

Bunch  Charge 

p  =  0.114  nC 

Bunch  Length  (FWHM) 

tb  =  0.4  ps 

Peak  Current 

I peak  =  285  A 

Normalized  Transverse  Emittance  (rms) 

£x,y  =  8  gm 

Energy  Spread  (rms) 

Ay/y=  0.4% 

Undulator: 

Period 

A  =  5.5  cm 

Length 

L  —  1.65  m 

Undulator  Parameter 

K  —  1.36 

Optical  Cavity: 

Wavelength 

A  =  1.6  gm 

Cavity  Length 

S  —  32  m 

Rayleigh  Length 

Zq  =  75  cm 

Quality  Factor 

Qn  —  5 

Mirror  Radius 

Rmir  —  3.8  Cm 

Table  5.1:  Jefferson  Laboratory  FEL  Oscillator  Parameters,  from  [17] 


d 

Figure  5.1:  Modeled  desynchronism  curve  (extraction  p  versus  dimensionless  desynchronism 
value  d )  for  Jefferson  Laboratory  14  kW  FEL  using  the  NPS  2D  model  (blue)  and  the  new  4D 
model  (red) 
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-42  x  42 

96%  in  cavity  fundamental  mode 

Figure  5.2:  Optical  wavefront  depiction  for  d  =  0.02.  Top:  optical  amplitude  |a|  versus  z  at 
x  =  y  =  0.  Bottom:  optical  amplitude  \a\  versus  (x,y)  at  z  =  0 


a  (0, 0, z)  15.3 


-57  x  57 

66%  in  cavity  fundamental  mode 

Figure  5.3:  Optical  wavefront  depiction  for  d  =  0.2.  Top:  optical  amplitude  \a\  versus  z  at 
x  =  y  =  0.  Bottom:  optical  amplitude  |a|  versus  (x,y)  at  z  =  0 
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The  2D  program  assumes  that  the  optical  wavefront  stays  in  the  fundamental  mode,  which 
leads  to  the  initial  conclusion  that  the  difference  between  the  curves  is  due  to  the  change 
in  optical  wavefront  mode.  Various  things  can  change  the  mode  of  an  optical  wavefront  in 
a  cavity,  including  mirror  imperfections  and  mirror  clipping.  Additionally,  gain  distortion 
can  cause  the  mode  to  change  [14].  The  2D  and  4D  model  runs  assumed  zero  mirror  imper¬ 
fections.  Also,  the  mirror  radius  used  for  the  modeling  was  sufficient  to  prevent  significant 
clipping  from  occurring.  One  difference  between  the  2D  and  4D  models  is  that  the  2D 
model  assumes  that  the  optical  wavefront  stays  in  the  fundamental  mode  and  the  4D  model 
allows  the  optical  wavefront  to  operate  in  higher-order  modes.  This  leads  to  the  hypothe¬ 
sis  that  the  mode  of  the  optical  wavefront  plays  a  role  in  the  width  of  the  desynchronism 
curve.  To  verify  this  conclusion  the  desynchronism  curve  was  generated  again,  except  the 
4D  parameters  were  adjusted  so  that  the  simulated  optical  cavity  prevented  the  higher  order 
modes  from  experiencing  significant  gain. 

Figure  5.4  shows  the  new  desynchronism  curve  where  the  2D  and  4D  curves  are  now 
aligned.  Investigating  the  optical  wavefront  at  d  =  0.02  and  d  —  0.2  shows  that  the  op¬ 
tical  wavefront  now  stays  in  the  fundamental  mode  (99%  and  97%,  respectively),  shown 
in  Figures  5.5  and  5.6.  These  results  support  the  conclusion  that  the  change  in  mode 
causes  the  model  desynchronism  curve  to  be  twice  as  long  as  a  curve  where  the  optical 
wavefront  stays  in  the  fundamental  mode. 
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Figure  5.4:  Adjusted  desynchronism  curve  for  Jefferson  Laboratory  14  kW  FEL.  For  the  4D 
model  (red),  the  cavity  mirrors  are  now  positioned  just  outside  the  undulatory,  forcing  the  optical 
wavefront  to  remain  primarily  in  the  fundamental  mode 
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a(x,y, 0) 
T=l.l 


-4x4 


99%  in  cavity  fundamental  mode 

Figure  5.5:  Adjusted  optical  wavefront  depiction  for  d  =  0.02.  Top:  optical  amplitude 
z  at  a:  =  y  =  0.  Bottom:  optical  amplitude  \a\  versus  (x,y)  at  z  =  0 


a  (0, 0, z)  15.4 


-4x4 


97%  in  cavity  fundamental  mode 

Figure  5.6:  Adjusted  optical  wavefront  depiction  for  d  =  0.2.  Top:  optical  amplitude 
z  at  x  =  y  =  0.  Bottom:  optical  amplitude  |a|  versus  (x,y)  at  z  =  0 
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Electron  Beam: 

Energy 

Eb  =  44.3  MeV 

Bunch  Charge 

q  —  5  pC 

Bunch  Length  (FWHM) 

tb  =  2.9  ps 

Peak  Current 

I peak  L85  A 

Normalized  Transverse  Emittance  (rms) 

£x,y  —  8.8  jU m 

Energy  Spread  (rms) 

Ay/y=  0.1% 

Undulator: 

Period 

A  =  3.3  cm 

Length 

L  =  4.95  m 

Undulator  Parameter 

K  =  0.709 

Optical  Cavity: 

Wavelength 

A  =  3.3  /i m 

Cavity  Length 

S  —  12.68  m 

Rayleigh  Length 

Zq  =  271.2  cm 

Quality  Lactor 

Qn  =  66.7 

Mirror  Radius 

Rjnir  ~  2  Cm 

Table  5.2:  Stanford  FEL  Oscillator  Parameters,  from  [19] 


5.2  Stanford  FEL 


The  Stanford  FEL  was  chosen  to  be  modeled  with  the  2D  and  4D  programs  since  it  was 
the  FEL  that  Benson  used  in  his  research  where  he  discovered  the  discrepancies  between 
experimental  data  and  modeling  results,  as  previously  discussed.  The  parameters  used  are 
shown  in  Table  5.2.  Figure  5.7  shows  the  desynchronism  curve  with  the  2D  program  in  blue 
and  the  4D  program  in  red  and  Figure  5.8  shows  the  experimental  desynchronism  curve. 
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1 


Figure  5.7:  Modeled  desynchronism  curve  for  Stanford  FEL 
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Figure  5.8:  Experimental  desynchronism  curve  for  Stanford  FEL,  from  [19] 


As  can  be  seen  in  Figure  5.7,  the  2D  and  4D  programs  match  and  give  an  approximate 
width  of  d  —  0.09  for  the  desynchronism  curve.  Analysis  of  the  optical  wavefronts,  as  seen 
in  Figures  5.9  and  5.10,  shows  that  it  stays  primarily  in  the  fundamental  mode  from  d  = 
0.0004  (100%  in  the  fundamental  mode)  to  d  =  0.06  (99%  in  the  fundamental  mode).  This 
is  further  verification  of  the  conclusion  made  in  the  last  section  regarding  diffraction  and 
the  correlation  between  the  optical  mode  and  the  width  of  the  desynchronism  curve.  The 
width  of  the  experimental  de synchronism  curve  in  Figure  5.8  corresponds  to  dimensionless 
desynchronism  d  ~  0.06,  according  to  (2.28).  This  width  is  shorter  than  the  width  of 
the  modeled  desynchronism  curve;  however,  it  is  on  the  order  of  half  the  width  of  the 
modeled  curve  reported  by  Benson  [19].  The  discrepancy  between  the  experimental  curve 
and  the  NPS  modeled  curve  may  be  due  to  the  fact  that  the  parameters  of  the  electron 
beam  were  slightly  different  than  that  reported  by  Benson;  it  may  also  indicate,  based  upon 
the  analysis  of  the  previous  section,  that  the  Stanford  experiment  had  higher  order  optical 
content  present  in  significant  amounts.  However,  the  NPS  code  does  come  significantly 
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closer  than  previous  modeling  in  predicting  the  width  of  the  desynchronism  curve. 


a(0,0,z)  156 


8x8 


1 00%  in  cavity  fundamental  mode 

Figure  5.9:  Optical  wavefront  depiction  for  d  =  0.004 


1 00%  in  cavity  fundamental  mode 

Figure  5.10:  Optical  wavefront  depiction  for  d  =  0.06 
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5.3  Emittance  Versus  Energy  Spread 

To  better  understand  the  robustness  of  the  design  of  a  notional  higher  power  FEL  oscillator, 
modeling  was  used  to  estimate  the  tolerances  of  the  emittance  and  energy  spread  of  the 
electron  beam  while  still  maintaining  a  decent  value  of  extraction.  The  parameters  of  the 
notional  oscillator  is  shown  in  Table  5.3. 

Figure  5.11  shows  a  3D  graph  of  extraction  ( 7] )  versus  transverse  emittance  (£x,y)  and  en¬ 
ergy  spread  (Ay/ y).  The  peak  extraction  achieved  is  2.7%;  however,  as  seen  in  Figure  5.11, 
this  value  is  only  achievable  at  an  emittance  of  5jum  and  an  energy  spread  of  0%  -  0.4%. 
Once  either  is  above  these  respective  levels,  the  extraction  has  a  significant  drop  off  and 
then  a  more  general,  and  predictable,  downward  trend. 


2.5 

2 

0.5 


Figure  5.11:  Extraction  (77)  for  emittance  ( eXJ )  versus  energy  spread  (Ay/y)  for  a  generic  oscil¬ 
lator 


Focusing  on  the  sharp  jump  in  extraction  shows  some  interesting  behavior  with  regard  to 
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Electron  Beam: 

Energy 

Eh  =  120  MeV 

Bunch  Charge 

q=  1  nC 

Bunch  Length  (FWHM) 

tb  —  2  ps 

Peak  Current 

1  peak  =  500  A 

Undulator: 

Period 

A  =  3  cm 

Length 

L  —  0.75  m 

Undulator  Parameter 

K=  1.8 

Optical  Cavity: 

Wavelength 

A  =  l.ljUm 

Cavity  Length 

S  =  20  m 

Rayleigh  Length 

Zo  =  7.5  cm 

Quality  Factor 

Qn  =  4 

Mirror  Radius 

Rmir  =  8  Cm 

Table  5.3:  Generic  FEL  Oscillator  parameters,  from  [17] 

the  optical  wavefront.  Figures  5.12  and  5.13  show  the  transverse  profile  of  the  electron 
pulse  and  optical  field  amplitude  in  the  x  direction  at  the  beginning  of  the  undulator  and 
the  evolution  of  the  optical  power  spectrum  at  the  pulse  center  for  eXJ  =  1 0  /i  m  and  5  /im, 
respectively,  with  Ay/y  =  0%  in  both  cases.  Examining  Figure  5.12  shows  that  the  FEL 
reaches  steady  state  with  an  extraction  of  1.7%.  The  optical  wavefront  is  a  Gaussian.  When 
the  emittance  is  changed  to  5/im,  the  FEL  finds  a  new  mode  of  operation  before  it  reaches 
steady  state.  Comparing  Figures  5.12  and  5.13  shows  this  new  transverse  mode.  The  optical 
field  amplitude  profile  has  developed  a  double  peak  and  the  optical  power  spectrum  shifts 
to  the  new  mode  prior  to  reaching  steady  state.  This  new  mode  dramatically  increases  the 
extraction  (now  2.7%).  This  new  mode  of  operation  is  not  seen  past  an  energy  spread  of 
0.4%  and  an  emittance  of  5/im  which  suggests  that  if  the  emittance  and  energy  spread  of  the 
electron  beam  can  be  maintained  below  these  limits,  then  extraction  will  see  a  significant 
increase.  While  these  results  suggest  a  behavior,  experiments  are  required  to  verify  the 
modeling  results  and  determine  the  implications  of  the  new  operating  mode. 

One  explanation  for  the  increase  in  extraction  is  that  for  the  lower  emittance,  the  electron 
pulse  has  a  higher  peak  current  at  the  center,  which  can  be  seen  when  comparing  Figures 
5.12  and  5.13.  This  higher  peak  current  will  increase  the  gain  of  the  optical  pulse  on 
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axis,  which  introduces  mode  distortion  by  placing  optical  power  into  higher  order  gaussian 
modes.  This  allows  the  cavity  to  store  additional  power  in  higher  order  modes  rather  than 
the  fundamental,  while  reducing  the  on  axis  optical  intensity  along  the  undulator,  perhaps 
delaying  the  onset  of  saturation.  While  the  behavior  provides  an  extra  percent  of  extraction, 
it  results  in  worse  beam  quality.  When  designing  an  FEL,  that  will  need  to  taken  into 
consideration. 
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Figure  5.13:  Optical  field  amplitude  and  optical  power  spectrum  for  £x.y  =  5/rm  and  A y/y  =  0% 


5.4  Mirror  Shift  and  Tilt 

To  better  understand  the  design  requirements  for  a  FEL  that  will  experience  shipboard 
conditions,  modeling  was  used  to  determine  the  effects  of  adjusting  the  mirror  shift  (A ym) 
and  mirror  tilt  (A0m)  on  extraction.  Previous  work  at  NPS  has  shown  that  an  FEL  will 
be  subjected  to  large  scale  but  low  frequency  vibrations  (due  to  waves  and  sea  state)  and 
smaller,  higher  frequency  vibrations  (due  to  onboard  machinery)  [20] .  The  larger  vibrations 
can  be  accounted  for  using  active  damping  that  is  commercially  available.  The  smaller,  high 
frequency  vibrations  are  more  difficult  to  mitigate  actively  and  can  become  a  problem.  The 
parameters  in  Table  5.3  were  used  with  an  emittance  of  10/im  and  an  energy  spread  of 
0.2%.  The  peak  extraction  is  1.7%  and  the  results,  as  seen  in  Figure  5.13,  show  that  at  least 
half  of  that  extraction  can  be  achieved  for  A ym  <  30/im  and  A0m  <  2/irad.  It  is  believed  by 
experimentalists  that  the  vibrations  that  are  experienced  during  shipboard  conditions  can 
be  actively  damped  down  to  support  the  mirror  shifts  and  tilts  to  maintain  half  of  the  peak 
extraction  for  the  modeled  FEL. 
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Figure  5.14:  Extraction  (rj)  for  mirror  shift  (A ym)  versus  mirror  tilt  (A 6m)  for  a  generic  oscillator 
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CHAPTER  6: 
Conclusion 


Analysis  of  desynchronism  curves  using  the  2D  and  4D  models  shows  that  the  model  curves 
are  longer  than  experimental  curves  when  the  optical  wavefront  stays  in  the  fundamental 
mode.  There  exist  a  few  possibilities  as  to  the  cause  of  this:  the  combined  effects  of 
diffraction  and  gain  distortion,  the  basis  set  for  the  optical  modes,  and  the  relative  lengths 
of  the  gain  medium  and  the  optical  cavity.  In  order  to  better  understand  this  issue,  more 
modeling  runs  would  be  required  as  well  as  experiments  that  simultaneously  measure  the 
desynchronism  curve  and  the  modal  composition. 

The  4D  simulation  results  used  to  determine  the  robustness  of  a  FEL  design  showed  that 
there  exists  a  new  transverse  mode  that  develops  for  a  normalized  emittance  £X;V  =  5  /im 
and  energy  spread  Ay/ y  <  0.4%.  Future  work  should  involve  the  determination  of  how 
this  new  mode  develops  as  well  as  experimental  verification.  If  this  new  mode  is  observed 
experimentally,  then  it  could  enable  a  significant  increase  in  the  extraction  and  will  help 
focus  the  efforts  of  FEL  design  to  achieve  the  required  electron  beam  parameters. 

Simulation  results  also  show  that  at  least  half  of  the  maximum  extraction  can  be  achieved 
with  a  mirror  shift  A ym  <  30/im  and  mirror  tilt  A0,„  <  2/irad.  These  results  show  the 
robustness  of  the  FEL  design  with  regard  to  mirror  alignment.  Further  work  should  focus 
on  expected  misalignments  due  to  vibrations  in  a  realistic  shipboard  environment.  Once 
that  is  determined,  active  mirror  alignment  systems  should  be  examined  to  improve  FEL 
performance. 

Modeling  continues  to  provide  valuable  insight  into  the  design  of  FELs  as  well  as  the 
dynamic  physics  of  FEL  operation.  Continued  modeling  will  be  required  to  further  advance 
FEL  technology  and  understanding  in  order  to  achieve  high  powers  for  operational  use. 
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